Adsorbed states of a long - flexible polymer chain 
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A phase diagram for a surface-interacting long flexible polymer chain in a two-dimensional poor 
solvent where the possibility of collapse exists is determined using exact enumeration method. A 
model of a self-attracting self avoiding walk (SASAW) on a square lattice was considered and up 
to 28 steps in series were evaluated. A new adsorbed state having the conformation of a surface 
attached globule is found. Four phases (i) desorbed expanded, (ii) desorbed collapsed, (iii) adsorbed 
expanded and (iv) surface attached globule are found to meet at a point on the adsorption line. 
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The behaviour of a long flexible polymer chain near 
an impenetrable surface is a subject of considerable ex- 
perimental and theoretical importance jl]-§| ■ This is be- 
cause such a chain may exhibit a phase diagram charac- 
terized by many different universality domains of critical 
behaviour. The subtle competition between the gain of 
internal energy and corresponding loss of entropy at the 
surface may lead to the possibility of the coexistence of 
different regimes and multicritical behaviour. An attract- 
ing surface may lead adsorption - desorption transition 
from the state when the chain is mostly attached to the 
surface, to the state of detachment when the tempera- 
ture is increased. This behaviour finds applications in 
lubrication, adhesion, surface protection etc. 

The essential physics associated with the behaviour of 
a surface interacting polymer chain in a good solvent 
where monomer- solvent attraction is greater than the 
monomer-monomer attraction, is derived from a model 
of self-avoiding walk (SAW) on a semi-infinite lattice. 
If the surface is attractive, it contributes an energy e a 
(< 0) for each step of the walk along the lattice bound- 
ary. This leads to an increased probability characterized 
by the Boltzmann factor u> = exp(— e a /kpT) of making a 
step along the wall, since for e a < 0, u> > 1 for any finite 
temperature T (kp is the Boltzmann constant). Because 
of this the polymer chain becomes adsorbed at low tem- 
peratures on the surface while at high temperatures all 
polymer conformations have almost same weight and non 
adsorbed (or desorbed) behaviour prevails. The transi- 
tion between these two regimes is marked by a critical 
adsorption temperature T a , with a desorbed phase for 
T > T a and adsorbed phase for T < T a . At T = T a 
one may define the crossover exponent <fi, as M ~ N^, 
where N is the total number of steps and M the num- 
ber of steps on the surface. The transition point T a is a 
tricritical point ||. Both the surface and the bulk crit- 
ical exponents have been calculated using renormaliza- 
tion group methods [Q , exact enumeration methods and 
Monte-Carlo simulations |^-||. For a two dimensional 
system exact values of the exponents have been found by 
using conformal invariance P]. 

The situation, is, however, different when the sur- 



face interacting polymer chain is in a poor solvent 
where monomer-monomer attraction dominates over the 
monomer- solvent attraction. As is well known, a long- 
flexible polymer chain in a poor solvent exhibits a tran- 
sition from a compact globule (collapsed state) to a ex- 
panded state when the temperature is increased. Above 
the critical 6 temperature (often referred to as the 9- 
point) the chain behaves as it would in a good solvent and 
below this temperature it behaves like a compact globule. 
At the #-point the chain behaviour is described by a tri- 
critical point of the 0(n) (n — > 0) spin system | jL0[ . How- 
ever, when the chain is in the vicinity of an impenetrable 
surface the competition between the monomer-monomer 
attraction and the surface-monomer interaction gives rise 
to many new features. Attempts have been made to 
study these features using several approaches P]7|,p ,pT| . 

For two-dimensions the transfer matrix method has 
been used for a directed polymer chain [|l2| whereas for 
the nondirected (isotropic) version the exact enumeration 
method has been used ||]. In both cases, three phases, 
desorbed expanded, desorbed collapsed and a single ad- 
sorbed phase have been predicted. However, the true 
nature of the phase diagram remained unknown. For 
three dimensions, the Monte-Carlo simulations method 
]l3| has been used for a finite length (~ 100) chain which 
led to a phase diagram containing four phases; desorbed 
expanded (DE), desorbed collapsed (DC), adsorbed ex- 
panded (AE) and adsorbed collapsed (AC). The phase di- 
agram shows a phase boundary between the AE and DC 
phases leading to two points on the phase diagram where 
three phases coexists (triple point). However, the phase 
diagram found by the exact enumeration technique JT^j] 
has many features which are different from that found in 
rcf. jl3). This indicates the possibility of a richer phase 
diagram than has been realized so far. In view of this we 
reexamine the problem of simultaneous adsorption and 
collapse of a linear polymer chain on a square lattice 
and investigate the phase diagram and critical param- 
eters using the exact enumeration technique. We prefer 
this technique because in this case the scaling corrections 
are correctly taken into account by a suitable extrapola- 
tion scheme. As shown by Grassberger and Hegger H, to 
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achieve the same accuracy with the Monte Carlo method 
one has to consider a polymer chain of about two or- 
ders of magnitude longer than in the exact enumeration 
method. 

We consider SASAW on a square lattice restricted to 
half space Z > (impenetrable hard wall). Walk starts 
from the middle of the surface. Let CN,N a ,N m be the 
number of SAWs with N steps, having N s (< N) step 
on the surface and N m nearest neighbor. We have ob- 
tained Cjv,iv s ,iv m f° r A < 28 for square lattice by exact 
enumeration method. 

Now we consider the interaction energy e a associated 
with each walk on the surface and e m for monomer- 
monomer interaction. Partition function of the attached 
chain is 
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where lo = e~ ta ^ kT and u = e~ tm l kT . lo > 1 and u > 1 
for attractive force. Reduced free energy for the chain 
can be written as 



G(co,u) = lim \ogZ N (uj,u) 

N— >oo iv 



(2) 



In general it is appropriate to assume that as N — » oo 

Z n {lo, u) ~ N^- V(w, u) N (3) 

where (J.(lo, u) is the effective coordination number and 7 
is the universal configurational exponents for walks with 
one end attached to the surface. The value of [i{lo, u) 
can be estimated using ratio method Jll| with associated 
Neville table. From equations (2) and (3) we can write 

log^(w,w) = lim — logZ N (uj,u) = G(u,u) (4) 

Zff{u, u) is calculated from the data of CN,N 3 ,N m using 
equation (1) for a given lo and u. From this we construct 
linear and quadratic extrapolants of the ratio of Z^{lo, u) 
for the adjacent values of N as well as the alternate one. 
Results for alternate N give better convergence. When 
u = 1 and lo = 1 the value of /i is found to be 2.638 
which is in very good agreement with the value given in 
ref. §,§. 

The value of lo c (u) at which polymer gets adsorbed for 
a given value of u is found from the (i) plot of G(lo,u) 
which remains fairly constants until lo = lo c and increases 
consistently as a function of lo, for lo > lo c (ii) from the 
plot of d 2 G(uj,u)/de a 2 at constant u and (m) from the 
plot of 7 — 71 ( see Eq.(6) and discussions which follow it) 
as a function of lo for different N. The value of lo c found 
from the plot of G(lo,u) is slightly lower than the peak 
value of d 2 G/de a 2 . It is, however, observed that as N is 
increased from 22 to 28 the peak value shifts to smaller 
lo and appears to converge on the value of lo found from 
G(oj, u) plot. We therefore choose the value of lo c found 
from the plot of G(oj,u) and determine lines ui c (u) and 



u>ci(u) (see Fig.l) by this method. For u = 1, the value of 
to c is 2.050 which is in very good agreement with the value 
(= 2.044 ±0.002) reported in ref. §]. Similarly the phase 
boundary separating the extended and collapsed phases 
is calculated from the plot of G(u>, u) as a function of u for 
a given lo. However, transition point u c is located more 
accurately from the peak of 8 2 G(lo, u)/de m 2 at constant 
lo. For u> = 1, the value of u c is 1.93 which is in good 
agreement with the value found by Foster et al [|| and the 
Monte Carlo results (= 1.94 ± 0.004) |J. The method is 
found to work for all values of lo i.e in both the bulk and 
the adsorbed regimes. However, as lo is increased the 
values of G(lo,u) do not remain as smooth as at lower 
values of lo, therefore introducing some inaccuracy in the 
value of u c . The estimate of this inaccuracy is of the 
order of 5% for lo > 4. We therefore conclude that the 
u c and lo C 2 (Fig.l) lines are determined with reasonable 
accuracy. 

The surface critical exponent 7^ ' k can be calculated 
using the relation: 
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where subscript "1" indicates the corresponding quantity 
of the surface (with one end of the polymer is attached to 
the surface). If we assume that [i does not depend on lo 
and polymer is in desorbed phase, then we can calculate 
the quantity 7 — 71 from the above equation 



7i = 



log{Z° N Z N _ 2 /Z N „ 2 Z N ) 



log(N/N 
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where superscript "0" indicates the corresponding quan- 
tity of the bulk (i.e without surface). In this case one 
calculates 7 — 71 for different N using above equation 
and plot it as a function of lo. The location of adsorp- 
tion point lo c can be determined from the intersection of 
successive approximation to 7 — 71 in the limit N — ► 00. 
This method, however, fails beyond the 0-point and re- 
produces closely the adsorption phase boundary below 
the 0-point found by the method discussed above. 

The phase diagram shown in Fig. 1 has four phase 
boundaries instead of three as reported in earlier work 
||. The u c line separates the expanded and collapsed 
phases. This line remains straight and parallel to w-axis 
in the bulk. This result is in agreement with that of ref. 
H . The special adsorption line lo c separates adsorbed 
expanded (AE) phase from that of desorbed expanded 
(DE). Beyond 0-point we have two boundaries w c \ and 
lo c2 . The line lo c \ separates the desorbed collapsed (DC) 
bulk phase from that of an adsorbed globule (AG) state, 
whereas the boundary lo c2 separates the AG phase from 
the AE phase. The point where u c line meets the special 
adsorption line lo c , all the four phases AE, DE, DC and 
AG coexist. The AG phase which exists between the 
boundaries lo c i and lo c2 for u > u c is essentially a two- 
dimensional globule sticking to the surface in the same 
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way as a liquid drop may lie on a surface. The existance 
of such a phase, to the best of our knowledge, is shown 
for the first time. 
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FIG.l: The phase diagram of a surface interacting linear 
polymer in 2-D space, ui and u axes represent, respectively 
the Boltzmann factor of surface interaction and monomer- 
monomer attraction. Regions marked by AE, AG, DE and 
DC represent, respectively, the adsorbed polymer in expanded 
( swollen) state and globular state, desorbed polymer in ex- 
panded and collapsed state. 

To find the monomer density as a function of distance 
normal to the surface for different regimes of the phase 
diagram we first found the terms which make most sig- 
nificant contribution in Eq.(l) for the walks of 26 steps. 
For these walks we evaluated the monomers (number of 
visited sites) on different lattice layers. The results are 
shown in Fig. 2 where we plot the fraction of monomers 
lying on different layers for many values of u and uj cor- 
responding to different parts of the phase diagram. Since 
the walks always start from the surface, there is at least 
one monomer on the surface for all the cases. In Fig. 
2(a) u is taken equal to 1.5 which is less than u c (1.93) 
and therefore it corresponds to expanded state. The val- 
ues of u) = 1.0 and 3.5 correspond, respectively to DE, 
and AE phases, while uj = 2.20 lies on the special ad- 
sorption line w c for u = 1.5. Fig. 2(b) shows the change 
in monomer density as uj is increased for u = 2.5. In 
this case uj = 1.0, 3.4 and 3.59 correspond to DC, AG 
and AE phases, respectively. For u = 2.5 uj c i = 2.75 and 
uj C 2 = 3.50, thus the value uj = 3.59 is just above the uj C 2 
line. The large change in monomer density distribution 
when uj value is changed from 3.4 (slightly below uj C 2 line) 
to 3.59 (just above uj C 2 line) is evident and confirms the 
existence of new (AG) phase. 

In Fig. 3 we plot the quantities < n 
giving the average fraction of monomers on the surface 
and number of pairs and defined as 
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FIG. 2: Fraction of monomers (n) on different layer. 

The transition points on each curve is marked by dot. 
While the dot on a curve of < n s > (in Fig. 3(a)) indicate 
the adsorption transition (i.e point on uj c and uj c1 lines 
depending on the values of u), the dot on a curve < n p > 
(in Fig. 3(b)) indicates the transition to collapsed state 
(i.e point on line u c and uj C 2 depending on the values of 
uj). We may note that the value of < n p > for AG phase 
is comparable to that in the DC phase. We have found 
that along the line uj c and oj c \, < n s >= 0.07 ±0.004 and 
along the line u c and uj c2 , < n p >= 0.5 ± 0.003. This 
shows that the error involved in determining the phase 
boundaries is very small. 

It is obvious from these results that when the chain gets 
adsorbed from the expanded bulk state (i.e for u < u c ), 
it acquires a conformation at uj ~ uj c (u) such that a small 
fraction (~ 10%) of monomers get attached to the surface 
(see Fig. 3(a)) and others are still in the bulk. Though 
the chain has formed a layer parallel to the surface there 
are considerable fluctuations in direction normal to the 
surface layer. As uj is increased for the same value of 
u, the fluctuations along the normal to the surface get 
suppressed and at large uj(» uj c ) the chain lies on the 
surface with very little fluctuations (see Fig. 2(a)). On 
the other hand, when the adsorbing chain was in col- 
lapsed bulk state, then at uj = uj c \ the collapsed chain 
gets attached to the surface. Here again the number of 
monomers getting attached to the surface are about 10% 
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(see Fig. 3(a)). For w c i < u> < uj C 2(u) the chain re- 
mains in the form of globule attached to the surface. In 
this range the monomer-monomer attraction remains ef- 
fective in holding the monomers in the neighbourhood of 
each other than the surface-monomer attraction whose 
tendency is to spread the chain on the surface (see Fig. 
2(b)). For lo > ujc2{u) the globule conformation becomes 
unstable as surface-monomer attraction becomes more ef- 
fective than the monomer-monomer attraction and there- 
fore the chain spreads over the surface (just like a liquid 
spreads over a wetting surface). 
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FIG. 3: (a) Average fraction of monomers on the surface 
(< n a >) as a function of u> and (b) number of pairs (< n p >) 
as a function of u. Dot on a curve indicates the transition 
point; in (a) adsorption - desorption and in (b) expanded - 
collapsed transitions. 

Summarizing we studied a SASAW in the presence of 
an attracting impenetrable wall and obtained the phase 
boundaries separating different phases of the polymer 
chain from data obtained by exact enumerations. We 
report a new adsorbed state which has conformation of 
a compact globule sticking to a surface in same way as 
a liquid drop may lie on a non wetting surface. The 
monomer density distribution, the number of monomers 
on the surface and the number of nearest neighbours in 
different regimes of the phase diagram are obtained. 

The work was supported by the Department of Science 
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